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Abstract 

A simple model for a granular ratchet corresponding to a single grain bouncing 
off a vertically vibrating sawtooth-like base is studied. Depending on the vibration 
strength, the sawtooth roughness and the restitution coefficient, horizontal trans- 
port in both the preferred and unfavoured directions is observed. A phase diagram 
indicating the regions in parameter space where each of the three possible regimes 
(no current, normal current, and current reversal) occurs is presented. 
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1 Introduction 



Ratchet-like motion, whereby directed transport is obtained from nonequilib- 
rium fluctuations, is important in many areas of scientific and technological 
interest, such as, molecular motors in biology and particle separation on nano- 
and microscales [1]. Recently, the problem of horizontal transport in a granu- 
lar layer vertically vibrated by an asymmetric base has also been considered, 
both experimentally [2] and computationally [3]. These so-called "granular 
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Fig. 1. A simple model for a granular ratchet. 

ratchets" have been shown to exhibit interesting dynamical features akin to 
those found in Brownian-motor models. 

In the present paper we approach the problem of granular ratchets by con- 
sidering the dynamics of a single grain bouncing off a vertically vibrating 
sawtooth-shaped base. The model proposed here is a generalization of a pre- 
vious model studied by some of us [4] for the gravity- driven motion of a par- 
ticle on an inclined rough surface with a staircase profile. In the present case, 
the 'rough surface' has a sawtooth shape with a horizontal baseline and is 
subjected to a vertical vibration. We find that, in spite of its simplicity, our 
single-particle model is able to reproduce some of the generic behaviors ob- 
served in experiments [2] and computer simulations [3] of granular ratchets, 
such as horizontal transport and current reversals. The main results of the 
paper are summarized in a phase diagram indicating the regions in parameter 
space where the different dynamical regimes (no current, normal current, and 
current reversal) can be observed. 



2 The model 



Our granular ratchet model consists of a particle bouncing off a vibrating 
sawtooth-like surface whose 'teeth' have a vertical height h and a horizontal 
length L, as shown in Fig. 1. A particle is launched with given initial velocity 
over the vibrating ratchet and then moves around through a succession of 
ballistic flights and inelastic collisions. We adopt a simple collision rule given 
by v' t = Vt and v' n = —ev n , where Vt and v n are respectively the velocity 
components tangential and normal to the collision plane, with prime denoting 
post-collisional velocities, and e is the coefficient of restitution. We assume 
furthermore that the sawtooth vibrates vertically with constant speed V = 
2A/t, where A is the vibration amplitude and r its period. For simplicity, we 
consider that A (and hence r) is negligibly small, so that the sawtooth base 
remains essentially at the same height, which allows us to determine explicitly 
the successive points of collision between the particle and the ratchet. The 
drawback here is that the effect of the vibration on the particle velocity is 
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taken into account in a probabilistic fashion, as will be discussed later. We 
believe, however, that the dynamics seen in our model is quite robust and 
should not be significantly altered if we consider a fully deterministic model. 

Suppose now that the particle takes off with velocity v = (u, v) from a point 
(x, y), which can be either on the sawtooth left edge (inclined ramp) or on its 
right edge (vertical wall). For convenience, we place our system of coordinates 
in such a way that the x axis is aligned with the inclined ramp and the y axis 
makes an angle = | — a, where a = tan _1 (L//i), with the vertical direction; 
see Fig. 1. (By convention, at each new flight we translate our system of 
coordinates to the tooth from which the particle departs.) Our goal then is 
to determine the next departure point (a/, y') and the corresponding takeoff 
velocity v 1 = (u',v'), so that the model dynamics can be described by a map 

F : (x,y,u,v) -> (x',y',u',v'). 



To this end, let us introduce the following dimensionless variables [4] 
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and the dimensionless roughness parameter k = s 2 /c 2 , where a = Lcosa, b = 
Lsina, c = cos0, and s = sin0. Hereafter we drop the capital notation with 
the understanding that we shall work solely with dimensionless quantities. 

Let us now define the jump number n G Z in such a way that its modulus 
equals the number of teeth the particle skips during the flight, with the con- 
vention that n is positive if the particle moves in the direction of the positive 
x axis. To calculate n, we must first determine the points (if any) where the 
particle trajectory crosses the line x + y — 1, connecting the tips of the saw- 
teeth (the dashed line in Fig. 1). Let be the y coordinates of these two 
crossing points (if they exist). An easy calculation shows [5] that 



± 1-x + Ky . ku + v± J(ku + v) 2 - (l- x -y)(l + K) 

Vh = ~TT^ 2k{u ~ v) 



If yf t £ (0, 1] (including the case where there is no real root) then the particle 
does not cross the line x + y — 1 and hence n — 0. On the other hand, if 
y^ G (0, 1] then the jump number is determined by the second crossing point 
y\. More precisely, n = —y\ , where \x~\ denotes the ceiling function (i.e., 
the smallest integer greater than x) . Having thus determined n, we then need 
to find the point at which the particle hits the surface. Note that if the particle 
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jumps to the right then it can land only on a ramp, but if the jump is to the 
left then the particle can hit either a ramp or a vertical wall. Let us first 
consider the case when the particle lands on a ramp. 

Collision with a ramp. Let (x*,y*) denote the collision point and (u*,v*) the 
particle velocity prior to the collision. Using the kinematics of ballistic motion 
and the fact that y* = —n, one obtains after a simple calculation that 

u* = u — v + \fv 2 + n + y, 

v* = -^Jv 1 + n + y, (1) 
x* = x — n(n + y) + 2k(u — v) (v + \fv 2 + n + y) . 



If the particle jumps to the left (this happens when u — v < 0), the condition 
n — x* < k must also be satisfied, otherwise the particle would have hit a 
vertical wall (this case will be treated below). On the assumption that the 
particle lands on the n-th ramp, the new departure point (x',y') and the new 
tangential velocity u' after the collision will then be 

x = x* — n, y' = 0, u' = u*. (2) 



Let us now calculate the normal velocity component v' after the collision. First 
note that if \v *| < V = cV (here we have dropped the subscript from v*), then 
the particle can never overtake the vibrating ramp and hence the collision will 
surely happen when the sawtooth is moving upwards. On the other hand, if 
\v*\ > V then the particle can hit the ramp either in its upward or downward 
motion. Furthermore, if the collision takes place during the downward motion, 
then depending on the particle velocity after the first collision there might 
(or not) be enough time for the sawtooth to revert its motion and hit the 
particle a second time. Under the assumption that the vibration amplitude 
is negligibly small, one can calculate the probabilities for all possible cases. 
Here we shall simply quote the final result and refer the interested reader to 
Ref. [6] for details of the calculation. If we define the parameters r = |f*|/U, 
Pi = (l + r _1 )/2, and P 2 = 1 — (l + e _1 )/r, then the outgoing velocity v' and 
the corresponding probability p are as follows: 

i) if re (0, 1] then v' = -ev* + (1 + e)V, p = 1; 

ii) if r £ (1,1 + e~ l ] then 

v > = -ev* + (1 + e)V, p = Pi, 
v 1 = eV + (1 + e) 2 U, p = l-P 1 ; 

iii) if r G (1 + e _1 , 1 + 2e _1 ) then 
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v > = -ev* + (1 + e)V, p = Pi, 
v> = - ev * - (1 + e)V, p = P 2 , 
v' = eV + (1 + e) 2 V, p=l-P 1 -P 2 ; 

iv) if r > 1 + 2c" 1 then 

v > = -ev* + (1 + e)V, p = Pi, 
v > = -ev* - (1 + e)V, p = l-P v 

This probabilistic collision rule together with (2) thus yield the map F for the 
case when the collision occurs with a ramp. Next we consider the case when 
the particle hits a vertical wall (i.e., u — v < and n — x* > k). 

Collision with a vertical wall. Here we first need to determine the actual col- 
lision point (xIvjUw) with the vertical wall and the particle velocity {u^.v^) 
just before such collision. It is easy to show [5] that 



x* w = x + 2nut* — Kt 



*2 



2vt* - 1* 2 , 



u 



u-t*, 



v-t*, 



(3) 



where 



_ K (n- 1) +n-x 

2k(u-v) ■ [ } 



According to our convention, the new departure point {x f ,y f ) is given by 



x x„ 



n. 



y =y w + n - 



(5) 
(6) 



Now, to determine the velocity (u', v') after the collision we simply apply the 
collision rule, v' t = v t and v' n = —ev n , since the vibration has no effect on 
the outgoing velocity in this case. To do that, we must first obtain (v t ,v n ) in 
terms of {u* w1 v^), apply the above collision rule, and then return to the (u', v') 
coordinates. Performing this calculation (details omitted), we find 



u' = (s 2 - ec 2 )u* w + c 2 (l + e)u* , 
v' = s 2 (l + e)u* w + (c 2 - es 2 )v* w . 



(7) 
(8) 



Equations (5)-(8) give the map F when the collision takes place with a vertical 
wall, thus completing the mathematical formulation of our model. 
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Fig. 2. Two-dimensional phase diagrams for the planes: (a) (V, k, e = 0.5) and (b) 
(V, k = 1, e). The regions with J = and J < are separated by a critical line, and 
the dark regions correspond to points where current reversal ( J > 0) occurs. 

3 Phase diagram 

Although the map F found above cannot be studied analytically owing to the 
presence of a stochastic term, its long-term dynamics can be easily investigated 
on the computer. To this end, let us define the current J in our model as the 
particle average horizontal velocity: J = (v x ), where the statistics is taken over 
a large number of collisions. In the present paper, we are mainly interested 
in mapping out the parameter space (V, k, e) according to which dynamical 
regime (J = 0, J < 0, or J > 0) is observed at a given point. For illustration 
purposes, it is best however to look at 2D slices of the parameter space and 
construct the corresponding phase diagram. Two such phase diagrams are 
shown in Figs. 2 for the planes (V, k, e — 0.5) and (V, k — 1, e). 

With the help of Fig. 2, the generic behavior of our ratchet model can be 
summarized, as follows. For given values of k and e, we see that for sufficiently 
small V there is no net current, i.e., J = 0, meaning that the particle remains 
confined in one of the "potential wells" of the sawtooth profile. Then as V 
increases past a critical value V c (K,e), we observe a transition to a regime 
with J < 0, where the particle, on average, moves horizontally in the direction 
favoured by the ratchet asymmetry, namely, to the left. In the two-dimensional 
phase diagrams of Fig. 2, the critical surface V c (k, e) appears as a critical line 
separating the region with J = from that with J < 0. Furthermore, for 
certain values of V, k and e, corresponding to the dark regions in Fig. 2, a 
horizontal flow in the 'unfavoured' direction (i.e., with J > 0) is established. 
Notice that these current-reversal regions appear as band-like structures, a 
fact that can be traced back to the periodic nature of the ratchet profile [5]. 
It is important to note, however, that current reversals can occur only for a 
sufficiently large k and sufficiently small e. This is clearly seen in Fig. 2, where 
there is no current reversal for k < 0.8 in Fig. 2a nor for e > 0.6 in Fig. 2b. 
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4 Conclusions 



We have studied a simple model for a granular ratchet corresponding to a 
single grain bouncing off a vertically vibrating sawtooth-like base. We observe 
that for small shaking velocity V, the particle remains confined between two 
teeth, with no net motion, whereas for larger V a nonzero horizontal current is 
established. In this case, for most values of the parameters the particle moves 
in the preferred direction, but current reversals also occur rather regularly. 
This complex behavior is summarized in the phase diagram presented in Fig. 2, 
which shows the regions in the parameter space (V, k, e) where each of the three 
possible regimes (no current, normal current, and current reversal) can occur. 
Our single-particle model is thus able to reproduce some of the characteristic 
behavior observed in experiments [2] and computer simulations [3] of granular 
ratchets, such as horizontal transport and current reversals. In particular, 
our model shows that for current reversal to occur two key ingredients are 
necessary: i) a sufficiently rough vibrating surface and ii) sufficiently large 
dissipation at collisions. Further investigation concerning the nature of the 
'phase transition' to the nonzero-current regime as well as a more detailed 
study of the current-reversal mechanism will be presented in a forthcoming 
publication. 
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